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Abstract. A chromospheric age distribution of 552 late- 
type dwarfs is transformed into a star formation history 
by the application of scale height corrections, stellar evolu- 
tionary corrections and volume corrections. We show that 
the disk of our Galaxy has experienced enhanced episodes 
of star formation at 0-1 Gyr, 2-5 Gyr and 7-9 Gyr ago, 
although the reality of the latter burst is still uncertain. 
The star sample birthsites are distributed over a very large 
range of distances because of orbital diffusion, and so give 
an estimate of the global star formation rate. These re- 
sults are compared with the metal-enrichment rate, given 
by the age-metallicity relation, with the expected epochs 
of close encounters between our Galaxy and the Magel- 
lanic Clouds, and with previous determinations of the star 
formation history. Simulations are used to examine the 
age-dependent smearing of the star formation history due 
to age uncertainties, and the broadening of the recovered 
features, as well as to measure the probability level that 
the history derived to be produced by statistical fluctua- 
tions of a constant star formation history. We show, with 
a significance level greater than 98%, that the Milky Way 
have not had a constant star formation history. 

Key words: stars: late-type - stars: statistics - Galaxy: 
evolution - solar neighbourhood 



1. Introduction 

The question whether the Milky Way disk has experienced 
a smooth and constant star formation history (hereafter 
SFH) or a bursty one has been the subject of a num- 
ber of studies since the initial suggestions by Scalo (1987) 



formation in the disk. The use of the word 'burst' for 
these features (usually lasting 1-3 Gyr) is based on the 
fact that all methods used to recover the SFH are likely 
to smear out the original data so that the star formation 
enhancement features could be narrower than they seem, 
or be composed by a succession of smaller bursts. In this 
s ense, they were named bursts A, B and C, after Majewski 
( |1993|) . 

The most efficient way to find the SFH is using the 
stellar age distribution, which can be transformed into a 
star formation history after various corrections. Twarog 
( 1980| ) summarized some of these steps. Although his SFH 
is usually quoted as an evidence for the constancy of the 
star formation in the disk, he states that during the most 
recent 4 Gyr, the SFH has been more or less constant, 
followed by a sharp increase from 4 to 8 Gyr ago, and a 
slow decline beyond that. His unsmoothed data were also 



reanalysed by Noh & Scalo (199C) who have found more 
signs of irregularity. 

Barry ( |1988[ ) has improved this situation substantially 
by using chromosphe ric age s. His conclusion was criticized 
by Soderblom et al. (1991), who showed that the empiri- 
cal data would be still consistent with a constant SFH if 
the chromospheric emission-age relation is suitably modi- 
fied. However, Rocha-Pinto & Maciel ( 1998| ) have r ecently 
argued that the scatter in Soderblom et al. (1991)'s Fig- 
ure 13, which is the main feature that could suggest a 
non-monotonic age calibration, is probably caused by con- 
tamination in the photometric indices due to the chromo- 
spheric activity. The chromospheric activity-age relation 



was also further investigated by Donahue (1993, 1998), 



and the new proposed calibration still predicts a non- 
constant SFH if applied to Barry's data. 
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3. In section 4, statistical significances for tlie SFH are pro- 
vided by means of a number of simulations. The impact of 
the age errors on the recovered SFH is also studied. Some 
comparisions with observational constraints are addressed 
in section 5, and each particular feature of the SFH is dis- 
cussed in section 6, in view of the results from the simu- 
lations and comparisons with other data. The case for a 
non-monotonic chromospheric activity-age relation is dis- 
cussed in section 7. Our final conclusions follow in section 
8. A summary of this work was presented in RPSMF. 

2. Converting age distribution into SFH 

Assuming that the sample under study is representative 
of the galactic disk, the star formation rate can be derived 
from its age distribution, since the number of stars in each 
age bin is supposed to be correlated with the number of 
stars initially born at that time. 

We use the same 552 stars with which we have de- 
rived the AMR (Paper I), after correcting the metallicities 
of the active stars for the mi deficiency (Gimenez et al. 



1991 ; Rocha-Pinto & Maciel 1998 ), which accounts for the 



influence of the chromospheric activity on the photomet- 
ric indices. The reader is referred to Paper I for details 
concerning the sample construction and the derivation of 
ages, from the chromospheric Ca H and K emission mea- 
surements. 

The transformation of the chromospheric age distribu- 
tion into history of the star formation rate comprises three 
intermediate corrections, namely the volume, evolutionary 
and scale height corrections. They are explained in what 
follows. 

2.1. Volume correction 

Since our sample is not volume-limited, there could be 
a bias in the relative number of stars in each age bin: 
stars with different chemical compositions have different 
magnitudes, thus the volume of space sampled varies from 
star to star. To correct for this effect, before counting the 
number of stars in each age bin, we have weighted each star 
(counting initially as 1) by the same factor used for 
the case of the AMR, where d is the maximum distance at 
which the star would still have apparent magnitude lower 
than a limit of about 8.3 mag (see Paper I for details). 

This correction proves to change significantly the age 
distribution as can be seen in Figure |^. 

2.2. Evolutionary corrections 

A correction due to stellar evolution is needed when a 




-with volume correction 
without volume correction 



Age (Gyr) 

Fig. 1. Chromospheric age distribution with and without 
volume correction, which was applied to our sample to 
allow the derivation of a magnitude-limited SFH. 




Mass (in solar masses) 

Fig. 2. Mass distribution of the sample. Masses were cal- 
culated from a mean mass-magnitude relation given by 
Scalo (198£). From the figure, we estimate a mass range 
of 0.8-1.4 Mq for our sample. Note the substantial ab- 
sence of massive stars, compared to the left wing of the 
mass distribution. The evolutionary corrections attempt 
to alleviate this bias. 



(Scalo pJ86| ). In Fi gure 3, the mass distribution is shown. 
We take the mass range of our sample as 0.8 to 1.4 Mq, 
which agrees well with the spectral-type range of the sam- 
ple from nearly F8 V to K1-K2 V. As an example for the 
necessity of these corrections, the stellar lifetime of a 1.2 
Mm is around 5.5 Gyr (see Figure below). This means 
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Mass (in solar masses) 



Fig. 3. Stellar Lifetimes from a variety of sources: Bahcall & Piran ( | 1983D ; VandenBerg ( |1985D , for Z = 0.0169 and 
Z = .0017; Eggleton, Fitchett & Tout ( |1989D ; Schaller et al. ( [1992[ ), for Z = 0.02 and Z = 0.001; Bressan et al. 
( |1993D , for Z = 0.02; Fagotto et al. (^99§i), for Z = 0.004. 



corresponds to i = 0), with mass between 0.8 and 1.4 Mq 
is 



N*{t) = i){t) \ (f>{m)dm, 



(1) 



where 0(m) is the initial mass function, assumed con- 
stant, and ipit) is the star formation rate in units of 
Mq Gyr^^pc^^. The number of these objects that have 
already died today is 



r(t) 



(f)(m) dm, 



(2) 



where mr{t) is the mass whose lifetime corresponds to t. 
From these equations, we can write that the number of 
still living stars, born at time i, as 



N°^''{t)^N*{t)-N\t). 



(3) 



Ar*(i) =e(t)iV°b>5(t), 



where 



e{t) 



1 



/5 



(5) 



(6) 



The number of objects initially born at each age bin 
can be calculated by using equation (^), so that we have to 
multiply the number of stars presently observed by the e 
factor. These corrections were independently developed by 



Tinsley (1974), in a different formalism. RPSMF present 



another way to express this correction in terms of the stel- 
lar lifetime probability function. We stress that all these 
formalisms yield identical results. 

The function rur (t) can be calculated by inverting stel- 
lar lifetimes relations. Figure ^ shows stellar lifetimes for 
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VandenBerg ( 1985| ) are underestimated probably due to 
the old opacity tables used by him. The agreement in the 
stellar lifetimes shows that the error introduced in the 
SFH due to the evolutionary corrections is not very large. 

The adopted turnoff-mass relation was ca lculated from 
the stellar lifetimes by Bressan et al. (1993) and Schaller 
et al. (1992), for solar metallicity stars: 



logrririt) = 7.59 - 1.25 logi + 0.05(logO^ 



(7) 



where t is in yr. This equation is only valid for the mass 
range SM© > m > 0.7 Mq. 

We have also considered the effects of the metallicity- 
dependent lifetimes on the turnoff mass. To account for 
this dependence, we have adopted the stellar lifetimes for 
different chemical compositions, as given by Bressan et al. 
( 1993|) and Fagotto et al. ( 1994 a, b). Equations similar to 
Eq. (|7|) were derived for each set of isochrones and the 
metallicity dependence of the coefficients was calculated. 
We arrive at the following equation: 



\ogmr{t) = a + blogt + c(logt)^. 



(8) 



where a = 7.62 - 1.56[Fe/H], b = -1.26 + 0.34[Fe/H], 
c = 0.05 - 0.02[Fe/H]. Since [Fe/H] depends on time we 
use a third-order polynomial fitted to the AMR derived 
in Paper I. In that work, we have also shown that the 
AMR is very affected at older ages, due to the errors in 
the chromospheric bins. The real AMR must be probably 
steeper, and the disk initial metallicity around —0.70 dex. 
The effect of this in the SFH is small. The use of a steeper 
AMR increases the turnoff mass at older ages, decreasing 
the stellar evolutionary correction factors (Equation 
As a result, the SFH features at young and intermediate 
age bins (ages lower than 8 Gyr) increases slightly related 
to the older features, in units of relative birthrate which 
is the kind of plot we will work in the next sections. 

Note that equation (||) does not reduce to equation 
(0) when [Fe/H] = 0. The former was calculated from an 
average between two solar-metallicity stellar evolutionary 
models, while the latter uses the results of the same model 
with varying composition. The difference in the turnoff 
mass from these equations amount 12-15% from 0.4 to 15 
Gyr. 

The initial mass function (IMF) also enters in the for- 
malism of the e factor. For the mass range under consid- 
eration, the IMF depends on the SFH, more specifically 
on the present star formation rate. It could be derived 
from open clusters, but they are probably severely affected 
by m ass segregation, unresolved binaries and so on (Scal o 
199^) . We have adopted the IMF by Miller & Scalo (|l979|) . 



- Z-dependent lifetimes 

- No Z-dependence, Miller-Scalo IMF 

- No Z-dependence. Salpeter IMF 



; (Gyr) 



Fig. 4. Stellar evolution correction factors. The curves 
stand for Equations (Q; dashed line) and (P; solid line) 
and Miller-Scalo's IMF. A third curve (dotted line) gives 
the results of using a Salpeter IMF with the turnoff-mass 
given by Equation (Q). 

(shown by dots) gives the results using a Salpeter IMF 
with the turnoff-mass given by Equation (^ . The e factor 
does not vary very much when we use a different IMF. Be- 
ing flatter than Salpeter IMF, the correction factors given 
by the Miller-Scalo IMF are higher. However, the effects 
of neglecting the metallicity-dependence of the stellar life- 
times are much more important in the calculation of this 
correcting factor. Since low-metallicity stars live less than 
their richer counterparts, the turnoff-masses at older ages 
are highly affected. In the following section, we will use the 
e factors calculated for metallicity-dependent lifetimes. 

2.3. Scale height correction 

Another depopulation mechanism, affecting samples lim- 
ited to the galactic plane, is the heating of the stellar 
orbits which increases the scale heights of the older ob- 
jects. To correct for this we use the following equations. 
Assuming that the scale heights in the disk are exponen- 
tial, the transformation of the observed age distribution, 
NQ{t), into the function N{t) giving the total number of 
stars born at time t is 



N{t) = 2H{t)Noit), 



(9) 



where H{t) is the average scale height as a function of 
the stellar age. A problem arises since scale heights are 
always given as a function of absolute magnitude or mass. 
To solve for this, we use an average stellar age correspond- 
ing to a given mass, f ollowi ng the iterative procedure out- 
lined in Noh & Scalo (|l990|) . This average age, (r), can be 
obtained by 



for a constant SFH, which gives an average value for the 



It) = 



J^-tN{t)dt^ 



(10) 
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Fig. 5. Scale heights given by Scalo ( 1986 , solid line) and 
Holmberg & Flynn ( ^000| , dotted line). 



average ages through the definition of H{t), equations 
and ® can only be solved by iteration. We use the chro- 
mospheric age distribution as the first guess No{t), and 
calculate the average ages (r)o- These are used to convert 
H{m) to H{t), and the star formation history is found by 
equation (|9|), giving Ni(t). This quantity is used to cal- 
culate (t)i and a new star formation rate, N2{t). Note 
that, in equation the quantity that varies in each it- 
eration is H{t), not the chromospheric age distribution 
NQ^t). Our calculations have shown that convergence is 
attained rapidly, generally after the second iteration. 

Great uncertainties are still present in the scale heights 
f or dis k stars. Few works have addressed them sinc e Scal o 
( |l986D 's review (see e.g., Haywood, Robin & Crezc |l997| ). 
We w ill be working with two different scale heights: Scalo 
( |1986D and Holmberg & Flynn ( ^000[ ), that are shown in 
Figure ||. Haywood et al.'s scale heights are just in the 
middle of these, so they set the limits on the effects in the 
derivation of the SFH. 

The major effect of the scale heights is to increase the 
contribution of the older stars in the SFH. Better scale 
heights would not change significantly the results, so that 
we limit our discussion to these two derivations. 



3. Star formation history in the galactic disk 

3.1. Previous chromospheric SFH determinations 

In Figure |^, we show a comparison between two SFHs, 
derived from chromospheric age distributions available in 



the l iterature: Barry ( |1988[ SFH given by Noh & Scalo 
1990| ) and Soderblom et al. ( |199l| , SFH given by Rana & 




Fig. 6. Comparison between ch romo spheric SFHs pub- 
lished in the literature: Barry ( |1988 , acco rding to Noh 
& Scalo 1990 ) and Soderblom et al. (1991, according to 
Rana & Basu 1992| ). T he po sition of bursts A, B and C 
(named after Majewski 1993) are marked. 



Note that the SFHs in Figure || are very similar to 
each other, a result not really surprising since Soderblom 
et al. have used the same sample used by Barry. On the 
other hand, the corresponding events in Barry's SFH ap- 
pears 1 Gyr earlier in Soderblom et al.'s SFH. The dif- 
ferent age calibrations used in these works are the sole 
cause of this discrepancy. Barry makes use of Barry et 



al. (1987)'s calibration which used a low-resolution index 



He 



r,+ fir 



analogous to Mount Wilson logi?^j^, while Soderblom et 
al. use a calibration derived by themselves. In F igure 0, 
we show a comparison of the ages for Barry (198S)'s stars 
using both age calibrations. The difference in the ages are 
clea rly ca used by the slopes of the calibrations. Barry et 
al. ( 1987 )'s calibration gives higher ages compared to the 
other calibration, which explains the differences in the cor- 
responding SFHs published. 

3.2. Determination of the SFH 

The three corrections described in section 2 are applied to 
our data in the following order: the age distribution is first 
weighted according to the volume corrections, then each 
age bin is multiplied by the e factor and we iterate the 
result according to equations (^) and (p^). The final re- 
sult is the best estimate of the star formation history. It is 
shown in Figure for an age bin of 0.4 Gyr and Scalo's 
scale height. There can be seen three regions where the 
stars are more concentrated: at 0-1 Gyr, 2-5 Gyr and 7-9 
Gyr ago. Beyond 10 Gyr of age, the SFH is very irregular, 
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4 6 8 10 

Soderblom et al. (1 991 )'s age 




Age (Gyr) 



Fig. 7. Comparison of stellar ages (in Gyr) in the calibra- 
tions by Barry et al. (1987)'s and Soderblom et al. (1991) 
The first age calibration seems to overestimate the chro- 
mospheric ages by around 1 Gyr. 



amplitude of the events. In this plot, the importance of 
the older bursts is increased, since in Holmberg & Flynn 
( 2000 ) the difference in the scale heights of the oldest to 
the youngest stars is greater than the corresponding value 
in Scalo's scale heights. 

We have used an extended nomenclature to that of 
Majewsky ( 1993| ) to refer to the features found. At the age 
range where bursts B and C were thought to occur double- 
peaked structures are now seen. Thus, we have used the 
terms Bl and B2, and CI and C2, to these substructures. 



Also shown is the supposed burst D, as Majewski (1993) 
had suggested. Their meaning will be discussed later. The 
lulls between the bursts were named AB gap, BC gap and 
so on. Some of us have previously referred to the most 
recent lull as 'Vaughan-Preston gap'. We now avoid the 
use of this term because: 



1. 



The Vaughan-Preston gap is a feature in the chromo- 
spheric activity distribution; 

Due to the metallicity-dependence of the age calibra- 
tion, the Vaughan-Preston gap is not linearly reflected 
in an age gap; 



3. Henry et al. (1996, hereafter HSDB) shows that the 
Vaughan-Preston gap is less pronounced than was ear- 
lier thought, and docs not resemble a gap but a tran- 
sition zone. 

Comparing with other studies in the literature, the 
SFH seems particularly different. There are still three ma- 
jor star formation episodes but their amplitude, extension 




6 8 10 12 

Age (Gyr) 



Fig. 8. Star formation rate for an age b in of 0.4 Gyr. The 

nomenclature used by Majewski ( 1993 ) was extended to 
be used with the main features of the SFH. The terms Bl 
and B2, and CI and C2, stand for substructures of the 
supposed bursts B and C, respectively. Also shown is the 
supposed burst D. The gaps between the peaks are named 
AB gap, BC gap, and so on. The upper and lower panels 



show the SFH using Scalo (1986) and Holmberg & Flynn 
( 2000| ) scale heights, respectively. 



with different scale heights. The first number refers to the 
SFH with Scalo's scale height, and the other refers to that 
with Holmberg & Flynn's. 

As we can see, the main events of our SFH seem to oc- 
cur earlier than the corresponding events in Barry's SFH, 
by approximately 1 Gyr. This can also be seen in Figure ^ 
the SFR from Soderblom et al. (1991)'s data have features 
earlier than Barry by about 1 Gyr. This comes mainly 



from the use of Soderblom et al. (1991)'s age calibration 
on which we have based our ages. This hypothesis is re- 
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This work 


Barry (1988) 


Number of 'bursts' 


3 


3 


Age of burst A 


0-1 Gyr 


0-1 Gyr 


Age of burst B 


2-5 Gyr 


4-6 Gyr 


Age of burst C 


7-9.5 Gyr 


7-11 Gyr 


Stronger burst 


B 


B 


Duration of the most recent lull (AB gap) 


1 Gyr 


< 3 Gyr 


(% of stars formed in burst A)/Gyr 


10.48/9.58 


8.92 


(% of stars formed in burst B)/Gyr 


10.40/9.72 


11.50 


(% of stars formed in burst C)/Gyr 


8.88/10.88 


11.92 



of occurrence comes from the shrinking of the chronologic 
scale of the age calibration. 

The narrowing of the AB gap is one of the main dif- 
ferences of our SFH and that found by Barry. This can be 
expected since our sample does not show a well-marked 
Vaughan-Preston gap , cont rary to what is fou nd in the 
survey of Soderblom ( 1985 ) , from which Barry ( |l988| ) se- 
lected his sample. 

Some other differences in the amplitude and duration 
of the bursts can be understood as resulting from the dif- 
ferences in the samples used by us and by Barry. Nearly 
70% of our stars come from HSDB survey. We h ave alr eady 
shown in Paper I that HSDB and Soderblom ( 1985 ) sur- 
veys have different chromospheric activity distributions. 
These are directly reflected in the SFH. 

We have found double peaks at bursts B and C. At the 
present moment we cannot distinguish these features from 
a real double-peaked burst (that is, two unresolved bursts) 
or a single smeared peak. However, it is interesting to see 
that the previous chromospheric SFHs give some evidence 
for a double burst C. In Figure || burst C also seems to be 
formed by two peaks. On the other hand, the same does 
not occur for burst B. The feature called B2 corresponds 
more closely to burst B in the previous studies, but at the 
age where we have found Bl, the other SFHs show a gap. 

The resulting SFH comes directly from the age distri- 
bution, in an approach which assumes that the most fre- 
quent ages of the stars indicate the epochs when the star 
formation was more intense. Both the evolutionary and 
the scale height corrections do not change the clumps of 
stars already present in the age distribution. The only cor- 
rection which could introduce spurious patterns in it is the 
volume correction, which must be applied before the other 
two. Figure |l] shows how it affects the age distribution. It 
is elucidating that the major patterns of the age distribu- 
tions are not much changed after this correction. We refer 




111, f,-. fV,C 



1 m-i 1 n rray 



fine 



1 n 



Age (Gyr) 

Fig. 9. Star formation rate calculated after disregarding 
outliers with weights exceeding 2a (solid line) and 3a 
(dashed line), after the volume corrections, compared to 
the history used throughout this paper. 



volume correction. In spite of it, it is necessary to know if 
the presence of stars with very high weights (due to their 
proximity and low temperature) could affect the results. 
Therefore, we have recalculated the SFH now disregard- 
ing the stars that have very high weights after the volume 
correction. We have cut the sample to those stars with 
weights not exceeding 2a and 3a. The resulting SFHs is 
compared to the SFH of the whole sample in Figure ^ It 
is possible to see that the presence of outliers does not af- 
fect the global result. The uncertainty introduced affects 
mainly the amplitude of the events, at a level similar to 
that introduced by the uncertainty in the scale heights. 
We believe that the volume correction has not impinged 
artificial patterns on the data, and that the star formation 
just derived reflects directly the observed distribution of 
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4-.1. Inconsistency of the data with a constant SFH 

There is a widespread myth on galactic evolutionary stud- 
ies about the near constancy of the SFH in the disk. This 
conies primarily from earlier studies setting constrai nts to 
the present relative birthrate (e.g., Miller & Scalo 1979 ; 
Scalo 1986 ). The observational constraints have favoured 
a value near unity, and that was interpreted as a constant 
SFH. 

This constraint refers only to the present star forma- 
tion rate. As pointed out by O'Connell (1997) and Rocha- 
Pinto & Maciel (1997), it is not the same as the star for- 
mation history. 

A typical criticism to a plot like that shown in Figure 
^ is that the results still do not rule out a constant SFH, 
since the oscilations of peaks and lulls around the unity 
can be understood as fluctuations of a SFH that was 'con- 
stant' in the mean. This is an usual mistake of those who 
are accustomed to the strong, short-lived bursts in other 
galaxies. 

The ability to find bursts of star formation depends on 
the resolution. Suppose a galaxy that has experienced only 
once a real strong star formating burst during its entire 
lifetime. The burst had an intensity of hundred times the 
average star formation in this galaxy, and has lasted 10^ 
yr, which are typical parameters of bursts in active galax- 
ies. Figure ^ shows how this burst would be noticed, in 
a plot similar to that we use, as a function of the bin size. 
In a bin size similar to that used throughout this paper 
(0.4 Gyr), the strong narrow burst would be seen as a fea- 
ture with a relative birthrate of 3.5. If we were to convolve 
it with the age errors, like those we used in Paper I, we 
could find a broad smeared peak similar to those in Fig- 
ure ||. For a biggest bin size (1 Gyr), the relative birthrate 
of the burst would be lower than 1.5. Hence, a relative 
birthrate of 1.5 in a SFH binned by 1 Gyr is by no means 
constant. A great bin size can just hide a real burst that, 
if occurring presently in other galaxies, would be accepted 
with no reserves. 

In the case of our galaxy, the bin size presently cannot 
be smaller than 0.4 Gyr. This is caused by the magnitude 
of the age errors. We are then limited to features whose 
relative birthrate will be barely greater than 3.0, espe- 
cially taking into consideration that the star formation 
in a spiral galaxy is more or less well distributed during 
its lifetime. Therefore, in a plot with bin size of 0.4 Gyr, 
relative birthrates of 2.0 are in fact big events of star for- 
mation. 

A conclusive way to avoid these mistakes is to calculate 
the expected fluctuations of a constant SFH in the plots 
we are using. We have calculated the Poisson deviations 

fnr a rnnst.p.nt .'-IPH rnmnnspd Viv .^.^9, stars Tn Picrnrp. 




Bin size (in Gyr) 

Fig. 10. The evanescence of a strong short-lived star for- 
mation burst due to the bin size. A star formation burst, 
lasting lO"^ yr, and with varying intensity (10, 50 and 100 
times more intense than the average star formation rate) 
was considered. The plot shows the value of the relative 
birthrate at the time of a burst It can be seen that for 
an age bin similar to that used throughout this paper, 
namely 0.4 Gyr, even the most strong and narrow burst 
would be represented by a feature not exceeding 3.5 in 
units of relative birthrate. 



smallest error bars correspond to Poisson errors (ivA^, 
where N is the number of stars in each metallicity bin). 
The thinner longer error bar superposed on the first shows 
the maximum expected error in the SFH, coming from 
the combination of counting errors, IMF errors and scale 
height errors. These last two errors were estimated from 
Figures]^ and^. The contribution of the scale height errors 
are greatest at an age of 3.0 Gyr, due to the steep increase 
of the scale heights around solar-mass stars. The effect of 
the IMF errors are the smallest, but grows in importance 
for the older age bins. 

From the comparison of the maximum expected fluc- 
tuations of a constant SFH and the errors in the Milky 
Way SFH, it is evident that some trends are not consis- 
tent with a constant history, particularly bursts A and B, 
and the AB gap. We can conclude that the irregularities 
of our SFH cannot be caused by statistical fluctuations. 

4-2. The uncertainty introduced by the age errors 

The age error affects more considerably the duration of the 
star formation events, since they tend to scatter the stars 
originally born in a burst. We can expect that this error 
could smear out peaks and fill in gaps in the age distribu- 
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Age (Gyr) 

Fig. 11. Star formation rate with counting errors. The 
error bars correspond to an error of ±\/N, where N is the 
number of stars found in each age bin. The dotted lines 
indicate the 2a variations around a constant SFR for a 
sample having 552 stars. The labels over the peaks are 
the same as in Figure |[ 

this age distribution is depopulated from old objects, since 
some have died or left the galactic plane. Our method to 
find the SFH makes use of corrections to take into account 
these effects. However, some features in the age distribu- 
tion could be caused rather by the incompleteness of the 
sample. These would propagate to the SFH giving rise to 
features that could be taken as real, when they are not. 

Thus, if we want to differentiate our SFH from a con- 
stant one, we must begin with age distributions, generated 
by a constant SFH, depopulated in the same way that 
the Galactic age distribution. With this approach, we can 
check if the SFH presented in Figure || can be produced by 
errors in the isochrone ages in conjunction with statistical 
fluctuations of an originally constant SFH. 

We have done a set with 6000 simulations to study this. 
Each simulation was composed by the following steps: 

1. A constant SFH composed by 3000 'stars' was built by 
randomly distributing the stars from to 16 Gyr with 
uniform probability. 

2. The stars are binned at 0.2 Gyr intervals. For each 
bin, we calculate the number of objects expected to 
have left the main sequence or the galactic plane. This 
corresponds to the number of objects which we have 
randomly eliminated from each age bin. The remaining 
stars (around 600-700 stars at each simulation) were 
put into an 'observed catalogue'. 

3. The real age of the stars in the 'observed catalogue' is 
shifted randomly according to the average errors pre- 



age of occurrence of the most prominent peak, ampli- 
tude and age of occurrence of the deepest lull. 

One of the problems that we have found is that due 
to the size of the sample, and the depopulation caused by 
stellar evolution and scale height effects, the SFH always 
presents large fluctuations beyond 10 Gyr. These fluctua- 
tions are by no means real. They arise from the fact that 
in the observed sample (for the case of the simulations, 
in the 'observed catalogue'), beyond 10 Gyr, the number 
of objects in the sample is very small, varying from to 
2 stars at most. In the method presented in the subsec- 
tions above, we multiply the number of stars present in 
the older age bins by some factors to find the number of 
stars originally born at that time. This multiplying fac- 
tor increases with age and could be as high as 12 for stars 
older than 10 Gyr; this way, by a simple statistical effect of 
small numbers, we can in our sample find age bins where 
no star was observed neighbouring bins where there are 
one or more stars. And, in the recovered SFH, this age 
bin will still present zero stars, but the neighbouring bins 
would have their original number of stars multiplied by 
a factor of 12. This introduces large fiuctuations at older 
age bins, so that all statistical parameters of the simulated 
SFHs were calculated only from ages to 10 Gyr. 

In Figure |l2|, we present two histograms with the sta- 
tistical parameters extracted from the simulations. The 
first panel shows the distribution of dispersions around 
the mean for the 6000 simulations. The arrow indicates 
the corresponding value for the Milky Way SFH. The dis- 
persion of the SFH of our Galaxy is located in the farthest 
tail of the dispersion distribution. The probability of find- 
ing a dispersion similar to that of the Milky Way is lower 
than 1.7%, according to the plot. In other words, we can 
say, with a significance level of 98.3%, that the Milky Way 
SFH is not consistent with a constant SFH. 

In panel b of Figure |l^, a similar histogram is pre- 
sented, now for the value of the most prominent peak that 
was found in each simulation. In the case of the Milky 
Way, we have Bl peak with 6 = 2.5. Just like the previous 
case, it is also located in the tail of the distribution. From 
the comparison with the values of the highest peaks that 
could be caused by errors in the recovering of an originally 
constant SFH, we can conclude with a significance level of 
99.5% that our Galaxy has not had a constant SFH. 



The use of Holmberg & Flynn (2000) scale heights in 
the simulations increases these significance levels to 100% 
and 99.9%, respectively. 

These significance levels refer to only one parameter of 
the SFH, namely the dispersion or the highest peak. For a 
rigorous estimate of the probability of finding a SFH like 
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Fig. 12. Distribution of parameters from 6000 simulations, using the scale heights from Scalo (1986). The left panel 
shows the dispersions around the mean SFH, while the right panel gives the value of the most prominent peak. In all 
the plots, the arrow indicates the corresponding value for the Milky Way SFH. 



box charts with the results of the 6000 simulations. Super- 
imposed on these box charts, we show the SFH, now cal- 
culated with Holmberg & Flynn (2000)'s scale heights. For 
the sake of consistency, the simulations shown in the figure 
also use these scale heights, but we stress that the same 
quantitative result is found using Scalo's scale heights. 

A lot of information can be drawn from this figure. 
First, it can be seen that a typical constant SFH would 
not be recovered as an exactly 'constant' function in this 
method. This is shown by the boxes with the error bars 
which delineate 2cr-analogous to those lines shown in Fig- 
ure |l^. The boxes distribute around unity, but shows a 
bump between 1 to 2 Gyr, where the average relative 
birthrate increases to 1.4. This is an artifact introduced by 
the age errors. In each individual simulation, the number 
of stars scattered off their real ages increases as a function 



older than it). This decreases the average star formation 
rate with respect to the original SFH, and the propor- 
tional number of young stars increases, because they are 
less scattered in age due to errors. This gives rise to a dis- 
tortion in the expected loci of constant SFHs. Note also 
the increase in the 2(T-region as we go towards older ages, 
reflecting the growing uncertainty of the chromospheric 
ages. 

The diagram allows a direct estimate of the probability 
for each feature found in the Milky Way SFH be produced 
by fluctuations of a constant SFH. The box charts gives 
the distribution of relative birthrates in each age bin. An 
average probability for the major events of our SFH are 
shown in Figure |l^, besides the features under interest. 
Rigorously speaking, the probability for the whole Milky 
Way SFH be constant, not bursty, can be estimated by the 
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Fig. 13. Box charts showing the results of the 6000 simulations, using Holmberg & Flynn ( ^OOOD 's scaleheights. The 
horizontal lines in the box give the 25th, 50th, and 75th percentile values. The error bars give the 5th and 95th 
percentile values. The two symbols below the 5th percentile error bar give the 0th and 1st percentile values. The two 
symbols above the 95th percentile error bar give the 99th and 100th percentiles. The square symbol in the box shows 
the mean of the data. Superimposed, the Milky Way SFH is shown. From the comparison with the distribution of 
results at each age bin, the probability to find each particular event in a constant SFH can be calculated. The numbers 
besides the major events give the probabilities for their being fluctuations of a constant SFH. 



of theses results show that the Milky Way SFH was by no 
means constant. 

4-.3. Flattening and Broadening of the Bursts 

Since the errors in the chromospheric ages are not negli- 
gible, a sort of smearing out must be present in the data. 
Due to this, a star formation burst found in the recovered 
SFH must have been originally much more pronounced. 
This mechanism probably affects much more older bursts, 
since the age errors are greater at older ages and the de- 
population by evolutionary and scaleheight effects is more 
dramatic. We can assume that if we found a feature like a 



do so, we use simulations of a SFH composed by a single 
burst over a constant star formation rate. The 'burst' is 
characterized by occurring at age r, having intensity c 
times the value of the constant star formation rate, and 
lasting 1 Gyr. We want to know the fraction of the burst 
that is recovered, as a function of age and of the burst 
intensity. 

We have performed 50 simulations for each pair (t, c) , 
with around 3000 stars in each simulation. A summary of 
these simulations is shown in Figure In all the panels 
(for varying c) , the fraction of the recovered burst is high 
for recent bursts and falls off smoothly until 8-9 Gyr, when 
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Fig. 14. Recovered fractions for a SFH composed of a single burst superimposed on a constant rate. The 'burst' is 
characterized by occuring at age r, having intensity c times the value of the constant star formation rate, and lasting 
1 Gyr. We show the cases for c = 1.5, 3.5, 5 and 10. In all the plots, the abscissa indicates the age t where the burst 
happened. The fraction recovered in the first 2 Gyr of age is greater than unity, due to the same problem that distorted 
the 2a loci of the constant SFHs in Figure ^ (see text). 



can conclude that it is more difficult to find bursts older 
than 8-9 Gyr, irrespective of its original amplitude. 

A second problem in the method is the broadening of 
the bursts. This depends sensitively on the age at which 
the burst occurs, and the results are even more dramatic. 
To illustrate this, another set of simulations was done. We 
consider now a SFH composed of a single burst, of 1000 
stars, lasting 0.4 Gyr. No star formation occurs except 
during the burst. We vary the age of occurrence from 0.3 
Gyr to 6 Gyr ago. Just one simulation was done for each 
age of occurrence, since we are only looking for the mag- 
nitude of the broadening introduced by the errors, so the 
exact shape of the recovered SFH does not matter. The 
recovered SFHs are shown in Figure ^ Only the younger 
bursts are reasonably recovered. The burst at 6 Gyr can 
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Fig. 16. Comparison between the metal-enrichment rate 
(defined as the derivate in time of the AMR, expressed 
by absolute metallicity Z), and the SFH. On theoreti- 
cal grounds both quantities should be correlated, and it 
should give a test to the reliability of the SFH. In practice, 
the magnitude of the errors in both functions still hinders 
the application of this test. 



5. Comparison with other constraints 

5.1. The SFH driving the chemical enrichment of the disk 

On theoretical grounds, there should be a correlation be- 
tween the SFH and the the age-metallicity relation (here- 
after AMR). The increase in the star formation leads to an 
increase in the rate at which new metals are produced and 
ejected into the interstellar medium. The correlation is not 
a one-to-one, since the presence of infall and radial flows 
can also affect the enrichment rate of the system. More- 
over, the enrichment rate is constrained by the amount 
of gas into which the new metals will be diluted. Never- 
theless, it is interesting to see whether the AMR we have 
found in Paper I is consistent with the SFH derived from 
the same sample, especially because, to our knowledge, 
this was never tried before. 

From the basic chemical evolution equations (Tinsley 



1980| ), for a closed box model (i.e., no infall), the link 



between the AMR and the SFH can be written as 

— it) cx — Vt, 

At^ ' m,g{t)' 



(11) 



where Z{t) gives the AMR, expressed by absolute metal- 
licity, ip(t) is the SFH as in equation (^, and mg{t) is the 
total gas mass of the system, in units of Mq pc^^. 

According to this equation, bursts in the SFH are 
echoed through an increase of the metal-enrichment rate. 



the other hand, the intrinsic metallicity dispersion of the 
interstellar medium can certainly somewhat obscure this 
proportionality, especially if it were as big as the AMR by 
Edvardsson et al. ( 1993| , hereafter Edv93) suggests. 

In Figure |l6|, we show a comparison between the metal- 
enrichment rate (top panel) with the SFH (bottom panel) . 
The enrichment rate increases substantially in the last 2 
Gyr, which could be a suggestion for a recent burst of SFH. 
However, the agreement between both functions seems 
very poor. There is a peculiar bump in the enrichment 
rate between 4 and 6 Gyr, which is coeval to a feature in 
the SFH, but most probably this is mere coincidence. 

Although we have used iron as a metallicity indica- 
tor, which invalidates Equation ([III), due to non recycling 
effects, we are not sure whether the situation would be 
improved by using O. The errors in both the AMR and 
SFH are still big enough to render such a comparison ex- 
tremely uncertain. However, it can be a test to be done 
with improved data. The more important result for chem- 
ical evolution studies is that, provided that we know ac- 
curately both functions, the empirical AMR and SFH will 
allow an estimate of the variation of the gas mass with 
time, which could lead to an estimate of the evolution of 
the infall rate. Future studies should attempt to explore 
this tool. 

5. 2. Scale length of the SFH 

The stars in our sample are all presently situated within a 
small volume of about 100 pc radius around the Sun. The 
star formation history derived from these stars is neverthe- 
less applicable to a quite wide section of the Galactic disk, 
since the stars which are presently in the Solar neighbour- 
hood have mostly arrived at their present positions from 
a torus in the disk concentric with the Solar circle. 

We have investigated how wide this section of disk is 
by integrating the equations of motion for 361 stars of 
the 'kinematic sample' (see Paper I) within a model of 
the Galactic potential. The potential consists of a thin 
exponential disk, a spherical bulge and a dark halo, and is 
described in detail in Flynn et al. ( 1996 ). For each star we 
determine the orbit by numerical integration, and measure 
the peri- and apogalactic distances, Rp and Ra and the 
mean Galactocentric radiu s, Rm = {Rp + Ra)/'2 for the 
orbit (cf. Edvardsson et al. 1993 ). 

The distribution of Rm is shown in Figure ^ Most of 
the stellar orbits have mean Galactocentric radii within 2 
kpc of the Sun (here taken to be at Rq = 8 kpc), i.e. 6 < 
Rm < 10 kpc. Very few stars in the sample are presently 
moving along orbits with mean radii beyond these limits. 



As discussed by Wielen, Fuchs and Dettbarn (1996) 
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5.3. The Galaxy and the Magellanic Clouds 




(kpc) 

Fig. 17. The histogram presents the mean galactocentric 
radius for the orbits of 361 stars presently located near 
the Sun. According to Binney & Sellwood (2000), these 
mean galactocentric radii indicate with good confidence 
the stellar birthplace radii. This shows that the SFH here 
derived is not a local history, but the average history over 
a significant part of the Milky Way disk. 



born i?,„(0). Widen, Fuchs and Dettbarn describe the pro- 
cess by which stars are scattered by these irregularities as 
orbital diffusion, and show that over time scales of several 
Gyr, that one cannot reconstruct from Rm the radius at 
which any particular star was born to better than a few 
kpc. This is of the same order as the width of the distribu- 
tion of Rm seen in Figure We therefore conclude that 
our stars fairly represent the star formation history within 
a few kpc of the present Solar radius, 6 < Rm < 10, or 
the "middle distance" regions of the Galactic disc. The 
SFH of the inncr-disk/bulge, and the outer disk are not 
sampled. 

However, Binney & Sellwood ( 200C| ) have criticized 
this conclusion. They show that during the lifetime of a 
star, the guiding-center of its orbit can change generally 
by no more than 5%. In this scenario, the value of Rm 
that we have calculated is close to the galactocentric ra- 
dius of the star birthplace, and our star formation history 
would still be representative of a considerable fraction of 
the galactic disk, 7 < Rm < 9. 

Another important conclusion of kinematic studies it 
that the older is a feature in the SFH, the more damped 
it is recovered from the data, related to its original ampli- 
tude (see, for example, Meusinger 1991bD , since the stars 
formed by the burst will be scattered through a larger re- 
gion. Hence, the younger bursts in our SFH are the most 
local features. This does not mean that they are most 
probably 'local irregularities'. In time scales of 1-2 Gyr, 



When evidences for an intermittent SFH in the Galaxy 
were first discovered, Scalo ( 1987 ) proposed that they 
could have originated from interactions between the 
Galaxy and the Magellanic Clouds. Indeed, the Magel- 
lanic Clouds are known to have probably experienced some 
episodes of strong star formation for a long time. Butcher 
( |1977| ) first proposed that the bulk of star formation in 
the Large Magellanic Cloud (LMC) has occurred from 3- 
5 Gyr ago, by the analysis of the luminosity function of 
field stars. Stryker et al. ( |l98lD and Stryker ( |1983| ) sub- 
sequently confirmed this result. In the last few years, ad- 
ditional studies have arrived almost at the same conclu- 



sions (Bertelli et al. |1992[ Vallenari et al. |l996a|b). West- 
erlund ( |199C ) also remarked that the star formation in 
the LMC seems to have been very small from 0.7 to 2 
Gyr ago. A very recent burst of star formation (around 
150 Myr a go) w as also found by the MACHO team (Al- 
cock et al. 1999) from the study of the period distribution 
of 1800 LMC ccpheids. Their analysis present compeling 
arguments favouring this hypothesis, as well as for the 
propagation of the star formation to neighbour regions. 

However, these results have more recently been ques- 
tioned, on the basis of colour-magnitude diagram synthe- 
sis. Some authors claim that important information on the 
SFH are provided by the part of the colour-magnitude 
diagram below the turnoff-mass, which could only be re- 
solved with the most recent observations (Holtzman et al. 
199£, and references therein; Olscn 19991 ). These pa- 



1997 



pers conclude that star formation in the LMC has been a 
continuous process over much of its lifetime. 

Note that continuity in the SFH does not means con- 



stancy. Holtzman et al. (1999) points that their method 



cannot constrain accurately the burstiness of the SFH in 
the LMC on small time scales, particularly for ages greater 
than 4 Gyr. Nevertheless, they show evidence for an in- 
crease in th e star formation rate in the last 2.5 Gyr. Dol- 
phin ( 2000| ) arrives to the same conclusion studying two 
different fields of the LMC, separated by around 2 kpc one 
from the other. The author recognizes that some large en- 
vironment alteration must have triggered an era of star 
formation in our neighbour galaxy. 

In spite of the controversy, it is impossible not to verify 
that some results on the SFH of the LMC are in appar- 
ently synchronism with some SFH events in the Milky 
Way disk. But this should be not really surprising. The 
Magellanic Clouds are satellites of our Galaxy, and past 
interactions between them were a rule, not an exception. 



Byrd & Howard (1992) showed that a companion satellite, 
whose mass is larger than 1% of the primary galaxy, could 
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Clouds can produce a dynamical wake in the halo that 
distorts the disk (Weinberg 1999[ ) . It is quite possible that 
such an effect could also enhance the star formation in the 
disk (M. Weinberg, private communication). 

Additional evidence comes from dynamical studies of 
the Magellanic Clouds. Several groups have worked on the 
derivation of their orbits around the Galaxy. The full or- 
bit of the Magellanic Clouds are still unknown, but there 
is some agreement in the published works. The most im- 
portant is that all of these works conclude that the most 
recent close encounter between the Clouds and the Milky 
Way has occurred 0.2-0.5 Gyr ago, which was the closest 
encounter through the entire history of the system (how- 
ever, Holtzman et al. 1999 mention an unpublished work 



by Zhao in which the last perigalacticon occurred 2.5 Gyr 



ago). Murai & Fujimoto (1980) calculated that other close 
encounters have occurred 1.5, 2.6 and 7.5 Gyr ago. Gar- 



diner et al. (1994) revisited Murai & Fujimoto (198C)'s 
model and recalculated the epochs of the close encounters 
as aro und 1 .6, 3.4, 5.5, 7.6 and 10 Gyr ago. However, Lin 
et al. (1995) have found different values: 2.6, 5.3, 8.4 and 
11.8 Gyr ago. 

From these results we can tentatively assume that, in 
the last 12 Gyr, the Clouds have had at most six close 
encounters with the Milky Way occurring more or less at 
0.2-0.5, 1.4-1.5, 2.6-3.4, 5.3-5.5, 7.5-8.4 and 10-11.8 Gyr 
ago. Some of these encounters are not predicted by all the 
authors, while some are in good agreement. For the sake 
of simplicity, we will refer to these encounters as I, II, III, 
IV, V and VI, respectively. 

There are similarities between the time of close en- 
counters and the events of our derived SFH. In Figure ^ 
we show the epoch of these encounters superimposed over 
our SFH. We can associate burst A with encounter I, peak 
Bl with encounter HI, and peak CI with encounter V. It 
is not unlikely that peak B2 could also be associated with 
encounter IV. On the other hand, encounter VI probably 
cannot be responsible for any of the features found beyond 
9 Gyr, since it occurs in an age range where the SFH is 
highly uncertain and subject to random fluctuations. 

A significant exception to the rule is encounter II. It 
is thought to have happened in the middle of the AB 
gap. It seems strange to think that a close encounter be- 
tween interacting galaxies could suppress the star forma- 
tion. Other mechanism should be responsible for the gap. 



On the other hand, Lin et al. (1995) have not found such 
an encounter. In fact, these authors predict that by this 
time, the Clouds would be located in their apogalacticon, 
more than 100 kpc away. 

Although the comparison is very premature, we con- 
clude that the data on the age distribution and orbits of 
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Fig. 18. Star formation history compared with the times 
of close encounters between the Galaxy and the Magellanic 
Clouds. 

properly answered, since there is still much uncertainty in 
the Magellanic Clouds close encounters, as well as on the 
chronologic scale of the chromospheric ages. 

6. The features of the Milky Way SFH 

We now can return to the discussion of the meaning of 
each feature found in the SFH derived in section 3. 

6.1. Burst A 

The most recent star formation burst is also the most 
likely burst to have occurred, since it has occurred in the 
very recent past, and so is less affected by the age errors. 
A recent enhancement in the SFH is also pres ent in nearly 
all pr evious investi gatio ns of the SFH (Sca lo 1987; Barry 
1988t Go mez et al. |1990|; No h & Scalo |199C| ; Soderbl om et 



al. 1991: Micela et al. 1993; Rocha-Pinto & Maciel 1997 



Chcreul et al. 1998| ), and is consistent with the distribution 



of sp ectral types in class V stars (Vereshchagin & Chupina 
|l993 ). It is not present in the iso chrone age distributions 
(Twarog 1980 Meusinger 1991a ) most probably due to 
the difficulty to measure ages for stars near the zero-age 
main sequence, where we expect to find the components 
of this burst in a HR diagram. 

We can conclude with confidence that it is a real fea- 
ture of the SFH. However, being the youngest, it is also 
the most local feature, because the younger stars have had 
no time to diffuse to larger distances from their birthsites. 
Thus, we cannot be sure (from out data only) whether 
this feature applies to the Milky Way as a whole. 

On the other hand, it is known that the Large Mag- 
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been closer than ever in their history (Lin et al. 1995| ). 
This suggests that burst A could be caused by tidal inter- 
actions between our Galaxy and the LMC. 



6.2. AB gap 

A substantial depression in the star formation rate 1-2 
Gyr ago was found by many studies, beginning with Barry 
( |1988| ; see also the SFH derived from the massive whit e 



1999| ). 



dwarf luminosity function derived by Isern et al. 
This gap appears, although not directly, in the chromo- 
spheric age distribution (the so-called Vaughan-Preston 
gap) and in the spectral type distribution, between A and 



F dwarfs (Vereshchagin & Chupina 1993). A quiescence 
between 1 and 2 Gyr is also visible in Chereul et al. ( 1998D , 
in their study of the kinematical properties of A and F 
stars in the solar neighbourhood. 

This feature has been present in all steps of our work, 
from the initial age distribution in Figure to the SFH. 
Note that the volume corrections have deepened this lull, 
but it has not changed its duration. 

The AB gap is likely to have lasted for a billion years. 
Previous studies have given a more extended duration for 
it, but we believe that it was caused by the use of a highly 
incomplete sample, together with a chromospheric age cal- 
ibration that does not account for the different chemical 
composition of the stars. Since it is a relatively recent fea- 
ture, it only samples birthsites over a radial length scale 
of 1-2 kpc. 

6.3. Burst B 

The small lull between the peaks Bl and B2 is not present 
in the initial age distribution (Figure appearing only 
after the volume corrections. It is very narrow, which could 
be most probably caused by hazardous small weights of 
the stars in these age bins, during the volume correction. 
This is why we have presently no means to distinguish 
burst B from a single burst or an unresolved double burst. 
At its age of occurrence, considerable broadening of the 
original features is expected. Either way, our simulations 
give strong support to this feature. 

Previous studies have found star formation enhance- 
ments around 4 Gyr ago (Scalo 1987; Barry 1988; 



Marsakov et a l. [19901 Noh & S calo |1990| ; Soderblom et al. 
1991| ; Twarog |1980| ; Meusinger |1991aD . Note that a strong 



concentration of stars around this age can also be found 
in the age distribution of Edv93's stars, that we show in 
Figure 

A significant exception is the SFH found by some of 
us (Rocha-Pinto & Maciel 1997). This paper suggests that 




Fig. 19. Age distribution of Edv93's sample. The concen- 
tration of stars around 4 Gyr shows evidence of burst B. 
The lack of objects with ages lower than 2 Gyr is not an 
evidence of the AB gap, since Edv93 avoided the inclu- 
sion of very young stars due to the difficulty to measure 
an isochrone age for them. 



several AMRs from the literature and different SFHs were 
found for each AMR . The SFHs recovered wit h the AMR 
from Twarog ( 198C| ) and Meusinger et al. (1991) were 



preferred compared to that found with Edvardsson et al. 



1993 (hereafter Edv93) AMR. To be consistent with our 



O TVl O I QT* 



present result, we need to compare the present SFH with 
that coming from Rocha-Pinto & Maciel's method for an 
AMR similar to that found from our sample (paper I). Our 
AMR now looks very similar to t he me an points of Edv93's 
AMR. Rocha-Pinto & Maciel ( |1997D have found, using 
Edv93's AMR, that Burst B could have around the same 
intensity as burst C, and also a narrow AB gap lasting 1 
Gyr at most. Figure [2C| shows a comparison between their 
SFH (for Edv93's AMR) and the present history binned 
by 1 Gyr intervals. 

6.^. BC gap and Burst C 

The existence of the BC gap is directly linked with how 
much credit we are going to give to Burst C. From Figure 
[Tsl , one could say that no burst could be found around 
8-9 Gyr, and all supposed features are artificial patterns 
created by statistical fluctuations. To reinforce this the- 
oretical expectation, we have done a simulation to show 
how the features above could be formed by a bursty SFH. 
We have considered a SFH composed by three bursts, one 
occurring at 0.3 Gyr, lasting 0.2 Gyr, and the other at 4 
Gyr, also lasting 0.2 Gyr, and the last ocurring at 9 Gyr, 
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Fig. 2 0. Comparison between Rocha-Pinto & Maciel 
( |l997D 's SFH (for Edv93's AMR) and the present 
chromospherically-based SFR binned by 1 Gyr intervals. 
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Fig. 21. Results from an example simulation to show how 
the observed features: Burst A and B and AB gap could be 
caused by two bursts separated by a lull of star formation. 
See text for description of the original SFH used in this 
simulation. 



The magnitude of the age errors prevents us from as- 
signing a good statistical confidence to this particular fea- 
ture. 

However, it is not implausible that we have overesti- 
mated the age errors. A decrease of 0.05 dex in the age 
errors could alleviate the situation and allow the identifi- 
cation of peaks (although highly broadened) younger than 
10 Gyr, which would suggest that burst C is a real feature. 
A better estimate of the age errors would not create new 
bursts, or flatten the recovered SFH in these age bins, 
but would give confidence limits for the ages where the 
features found are likely to be real and not just artifacts. 

6.5. Burst D 



The so-called burst D was proposed by Majewski (1993) 



as a star formation event that would be responsible for 
the first stars of the disk, before the formation of the thin 
disk. 

A superficial look at Figure ^ could give us the impres- 
sion that the peaks beyond 11 Gyr were remnants of this 
predicted burst. However, as we have shown above, it is 
presently impossible to recover the SFH correctly at this 
age range, even if our age errors are overestimated by as 
much as 0.05 dex. The SFH at older ages are dominated by 
fluctuations, superimposed on the original strongly broad- 
ened structures, in such a way that it is imposible to dis- 
entangle statistical fluctuations from real star formation 
events. 

Theoretically, patterns as old as 13 Gyr could be found 
in the SFH, provided that they occurred not very close to 
younger ones, if the age errors were decreased by 0.10 dex, 
but that is hardly possible to be attained at the present 
moment since it would need to be of the order of magni- 
tude of the error in the logi?^j^ index. 

For these reasons, we give no credit to the peaks be- 
yond 11 Gyr in Figure ^. If burst D has ever occurred, 
probably the present chromospheric age distribution is not 
an efficient tool to flnd its traces. 



in Figure ^ Although the two more recent bursts can be 
well recovered, there is no sign of burst C at 9 Gyr. We 
have tried other combinations between the amplitude and 
time of occurrence of them, but in all cases the stars of 
burst C were much scattered from its original age. 

If on theoretical grounds there is no convincing argu- 
ments to accept the existence of burst C, the same does not 
occur on observational grounds. This puzzling situation 
comes from the fact that burst C has appeared in a num- 



7. The shape of the chromospheric activity age 
relation 

Soderblom et al. ( |l99lD argued that the interpretation of 
the chromospheric activity distribution as evidence for a 
non-constant SFH is premature. Particularly, the authors 
have shown that the observations do not rule out a non- 
monotonic chromospheric activity-age relation, even con- 
sidering that the simplest flt to the data is a power-law, 
like the one we used. 

Presently, there is good indication that the chromo- 
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Fig. 22. Comparison between the chromospheric activity- 
age relation we have used throughout this paper (solid 
line), the constant-sfr calibration proposed by Soderblom 
et al. (dotted line) and the constant-sfr calibration from 
our sample (dashed line). T he da ta and symbols are the 
same from Soderblom et al. ( 1991 ): 0, the sun; H, Hyades; 
U, Ursa Major Group; open diamonds, field F dwarfs; 
filled diamonds, binaries. 



relation, but the scatter is such that it is not presently 
possible to know whether the chromospheric activity de- 
creases steadily with time, or there are plateaux around 
some 'preferred' activity levels. There is a possibility that 
the clumps we are seeing in the chromospheric age distri- 
bution (which are further identified as bursts) are artifacts 
produced by a monotonic chromospheric activity-age re- 
lation. 

To keep the constancy of the SFH, Soderblom et al. 
( |l99lD proposed an alternative chromospheric activity- 
age relation that is highly non-monotonic. We have 
checked this constant-sfr calibration with our sample, but 
the result is not a constant sfr. This is expected, since 
there are many differences in the chromospher ic sam ples 
used by Soderblom (|l985| ) and Soderblom et al. ( 1991 ) and 
the one we have used (see our Figure 11 in Paper I). We 
have calculated a new constant-sfr calibration, in the way 
outhned by Soderblom et al. (1991). We have used 328 
stars from our sample (just the stars with solar metallic- 
ity, to avoid the metallicity dependence of logi?^^:), with 



ibration proposed by Soderblom et al. (dotted line) and 
the constant-sfr calibration from our sample (dashed line). 
The data and symbols are the same from Soderblom et 
al. ( |1991| ) . Both constant sfr calibrations agree reasonably 
well for the active stars, but deviate somewhat for the in- 
active stars. This is caused by the fact that to be consistent 
with a constant sfr, the calibration must account for the 
increase in the relative proportions of inactive to active 
stars, especially around logi^^K = —4.90, after the sur- 
vey of HSDB. Note that, our constant-sfr chromospheric 
activity-age relation is still barely consistent with the data 
and cannot be ruled out. There are few data for stars older 
than the Sun in the plot, and it is not possible to know 
whether the plateau for logi?^K < —5.0 in this calibration 
is real or not. 

We acknowledge that, given no other information, it 
is a subjective matter whether to prefer a complex star 
formation history or a complex activity-age relation. Nev- 
ertheless, there are numerous independent lines of evi- 
dence that also point to a bursty star formation history; 
the most recent and convincin g is th e paper by Hernan- 
dez, Valls-Gabaud & Gilmore ( ^OOOl) . They use a totally 
different technique (colour-magnitude diagram inversion) 
and find clear signs of irregularity in the star formation. 
In section 6, we listed several other works that indicate 
a non-constant star formation history, and the majority 
of them use different assumptions and samples. Strongly 
discontinuous star formation histories are also found for 



some galaxies in the Local Group (see O'Connell |1997D , 
in spite of the initial expectations during the early studies 
of galactic evolution that these galaxies should have had 
smooth star formation histories. 

For all these methods to give the same sort of re- 
sult, all the different kinds of calibrations would have to 
contain complex structure. It is simpler to infer that the 
star formation history is the one that actually has a com- 
plex structure. We think that when several independent 
methods all give a similarly bursty star formation history 
(although with different age calibrations, so they do not 
match exactly), our conclusion is supported over the irreg- 
ular activity-age but constant star formation rate solution. 

8. Conclusions 

A sample composed of 552 stars with chromospheric ages 
and photometric metallicities was used in the derivation 
of the star formation history in the solar neighbourhood. 
Our main conclusions can be summarized as follow: 

1. Evidence for at least three epochs of enhanced star 
formation in the Galaxy were found, at 0-1, 2-5 and 
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parameters, although the use of Fe as a metaUicity in- 
dicator, and the magnitude of the errors in both func- 
tions can still hinder the test. 

3. We examined in some detail the possibility that the 
Galactic bursts are coeval with features in the star 
formation history of the Magellanic Clouds and close 
encounters between them and our Galaxy. While the 
comparison is still uncertain, it points to interesting 
coincidences that merit further investigation. 

4. A number of simulations was done to measure the 
probability for the features found to be consistent with 
a constant SFH, in face of the age errors that smear 
out the original features. This probability is shown to 
decrease for the younger features (being nearly 0% for 
the quiescence in the SFH between 1-2 Gyr) , such that 
we cannot give a strong assertion about the burst at 
7-9 Gyr. On the other hand, the simulations allow us 
to conclude, with more than 98% of confidence, that 
the SFH of our Galaxy was not constant. 

There is plenty of room for improvement in the use of 
chromospheric ages to find evolutionary constraints. For 
instance, a reconsideration of the age calibration and a 
better estimate of the metallicity corrections could dimin- 
ish substantially the age errors, which would not only im- 
prove the age determination but also give more confidence 
in the older features in the recovered SFH. 
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